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We analyze the electrical characteristics of a circuit consisting of a free thin-film magnetic layer and 
source and drain electrodes that have opposite magnetization orientations along the free magnet's 
two hard directions. We find that when the circuit's current exceeds a critical value there is a sudden 
resistance increase which can be large in relative terms if the currents to source or drain are strongly 
spin polarized and the free magnet is thin. This behavior can be partly understood in terms of a 
close analogy between the magnetic circuit and a Josephson junction. 
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I. INTRODUCTION 

Electronic transport can usually be described in terms 
of effectively independent electrons. Recently, with 
the discovery and exploitation of spin-transfer torque^— 
(STT) effects, magnetism has joined superconductivity 
as an instance in which collective and quasiparticle con- 
tributions to transport are entwined in an essential way. 
The similarity between the non-equilibrium properties 
of magnetic and superconducting- systems is especially 
close when comparing the properties of a superconduct- 
ing circuit containing a Josephson junction to a magnetic 
circuit containing a free ferromagnetic layer with strong 
easy-plane anisotropy. As we explain in more detail be- 
low, the role of the Josephson junction bias current in the 
superconducting circuit is played in the magnetic case 
by the component of the spin-current injected into the 
nanoparticlc that is perpendicular to the easy plane. 

The electrical properties of a circuit containing a 
Josephson junction typically change drastically when the 
junction's critical current is exceeded. In this paper we 
propose that the magnetic circuit illustrated in Fig. [TJ 
which we refer to as a magnetic Josephson junction 
(MJJ), can exhibit similar drastic effects when a critical 
current related to the free magnet's in-plane anisotropy is 
exceeded. We show that the resistance of an MJJ can in- 
crease substantially when its critical current is exceeded 
provided that either source or drain currents are strongly 
spin polarized and magnetization relaxation in the free 
magnet layer is not too rapid. 

The analogy we explore is closely related to the early 
suggestion by Berger— that 180° domain walls in ferro- 
magnets should exhibit behavior analogous to the be- 
havior of Josephson junctions. Indeed the MJJ geome- 
try we propose may be viewed as a discrete version of 
a pinned 180° domain wall. Although the magnetiza- 
tion dynamics induced emf that Berger predicted based 
on an analogy to the ac Josephson relation has now^ 
been confirmed experimentally, electrical signals of mag- 
netization texture dynamics in uniform bulk materials 
tend to be weak. The MJJ geometry we study is also 
closely related to that employed in spin-transfer torque 
oscillators.- It is well known that the dc resistance of an 
STT-oscillator tends to increase once the magnetization 
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FIG. 1: (color online). A magnetic Josephson junction con- 
sists of a free magnetic layer with a well-defined easy-plane 
and weak in-plane anisotropy separated by non-magnetic 
spacers from magnetic source and drain electrodes whose 
magnetizations are opposite and perpendicular to the easy 
plane. In this article we choose the spin z-axis along a free 
magnet hard direction. 

becomes dynamic^ The increase in bias voltage at fixed 
current is closely related to Berger's Josephson voltage. 
From this point of view, the issue we explore theoret- 
ically in this paper is the possibility of achieving large 
relative changes of resistance in an STT-oscillator when 
the magnetization becomes dynamic. We believe that 
such an examination is motivated by recent advances in 
ferromagnetic metal spintronic o 10 ' 11 which have enabled 
the injection of more highly spin polarized currents and 
decreased the thickness of free magnetic layers, making 
them easier to manipulate electrically. MJJ devices not 
only provide a test bed for our understanding of the fun- 
damental physics of spin gauge fields, but could also be 
useful because of their unusual transport properties. 

Our paper is organized as following. In Sec. II, we 
comment at greater length on the relationship between 
Josephson junctions and easy-plane free magnets. In Sec. 
Ill, we discuss the theoretical framework we used for an- 
alyzing the transport properties of an MJJ. In Sec. IV 
and Sec. V, we identify two regimes in which the MJJ 
could have quite different resistances. Finally, in Sec. VI, 
we summarize our results and discuss our conclusions. 



II. NANOMAGNET EQUATIONS OF MOTION 

As shown in Fig.[IJ the MJJ is a multilayer heterostruc- 
ture consisting of ferromagnetic and non-magnetic nor- 
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mal metal layers. The two ferromagnetic electrodes have 
opposite magnetization directions, serving as spin polar- 
ized current source and drain. The free magnetic layer in 
the middle is sandwiched between two non-magnetic nor- 
mal metal spacer layers to avoid direct magnetic coupling 
with the electrodes. We assume that the free magnetic 
layer is small enough that its magnetization is spatially 
constant. Its magnetization direction dynamics is then 
described by a Landau-Liftshitz-Gilbert (LLG) equation, 



equivalent to the resistively and capacitively shunted 
Josephson junction models with easy-plane anisotropy 
playing the role of the capacitance, 
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in-plane anisotropy playing the role of the superconduct- 
ing coupling which supports a supercurrent 
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where m is a unit vector along the magnetization direc- 
tion of the free layer, 7 = \g\^B/h > is (minus) the gy- 
romagnetic ratio, H e g is an effective magnetic field, a is 
the dimensionless Gilbert damping constant and the right 
hand size of the equation above denotes a Slonczewskiir— 
spin-transfer torque termi^ that captures the coupling 
between transport and collective dynamics: 



dm 
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(rh x I 8 ) x rh. 
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Here S = M s V/j is the magnitude of the total spin of the 
free magnetic layer, V is the free layer volume, M s is its 
magnetization per unit volume, I s is the net spin current 
flowing out of the free layer and (rh x I s ) x rh = I s — 
(I s ■ m)m selects the component of I s that is transverse 
to rh. We assume that there is no applied magnetic field. 
The effective magnetic field, 
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then arises mainly from the magnetostatic energy of a 
thin film with an elliptical cross-section: 



K 7 ml + Kr, sin 2 
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where <j> is the magnetization orientation azimuthal an- 
gle and the anisotropy parameter K z 3> K p so that the 
magnetization direction is always close to its easy plane. 
When the magnetization direction is expressed in terms 
of m z and 0, and m z is assumed to be small, the LLG 
equations take the form 
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As discussed later we will be most interested in circum- 
stances under which rh z = and if, the component of 
spin-current in the azimuthal direction, vanishes. In that 
limit the right hand side of the equation for <j> reduces to 
the difference between the chemical potential of majority 
and minority spin electrons. Eq. ([5]) is then precisely 



'sup 



K p sin(20), 



(8) 



atj) playing the role of a dissipative quasiparticle current, 
and the spin-transfer torque term playing the role of junc- 
tion bias current. 

It is well known that the electrical properties of a 
Josephson junction changes drastically at its critical cur- 
rent. Due to the connections revealed above, we expect 
that similar behavior can occur in an MJJ. Because the 
MJJ bias current is a spin-current not a charge current, 
and its dissipationless supercurrent transports spins from 
up to down states not charges from one side of the junc- 
tion to the other, the current-voltage relationship of a 
circuit in which an MJJ is embedded is more subtle than 
in the superconducting Josephson junction case, and col- 
lective effects are generally weaker. In the following, we 
analyze MJJ I-V relationships with the aim of identifying 
circumstances under which the resistance change at the 
critical current can be large in circuits that are fabricated 
with modern spintronic materials. 



III. MAGNETOELECTRONIC CIRCUIT 
THEORY 

For small bias voltages the nanoparticle magnetization 
can achieve a static value by balancing the circuit bias 
and collective transport terms: l z s — I SU p- Once J sup 
reaches its maximum value, / cr ;t = K p , the magnetiza- 
tion's azimuthal angle becomes time-dependent and the 
time average of I sup quickly becomes negligible. In this 
dynamic limit must be balanced against the dissipa- 
tive term — a<p. Our main goal in this section is to shed 
light on the requirements for a substantial relative change 
in circuit resistance when the MJJ's critical current is ex- 
ceeded. For this purpose we use magnetoelectronic cir- 
cuit theory— _ — which simplifies the transport problem 
in ferromagnet/normal metal heterostructures by drop- 
ping some non-local effects and separating the circuit into 
normal metal and magnetic constituents within which 
the exchange field is taken to be constant in direction. 
In the following, we assume that the ferromagnetic layer 
is thicker than the magnetic coherence length^ which is 
on the order of A in transition metals. 

Magnetoelectronic circuit theory's key equation relates 
the charge and spin currents that flow from a ferromag- 
net into a normal metal to the ferromagnet's chemical 
potential and to the chemical potential of the normal 
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metal, including its non-equilibrium spin-accumulation is a dimensionless spin-dependent conductance, /xjvc is 
contributions. The charge current the normal metal chemical potential averaged over spin 

states while [iNs represents the magnitude and direction 
I c = — [2(3 1 "' 1 ' + g^) (unc — Mf) + (.9^ - 9^) /*jv s ■ riiFjh 1 spin-space of the spin-dependent part of the normal 
^ /g\ metal's chemical potential. For static magnetization the 

• ■ ... 3jI13j1o£OX1S ST31H - CUrrGIlt GXT)I"GSS10r^^™ IS 

Here f and 4- refer to majority and minority spins^ e> i J y 

m F is the ferromagnet's magnetization direction, g a,<7 



I s = -— [2(g t ' t - g^ifiNc ~ MF)m F + (,g t ' t + g i - l )((J, Ns ■ m F )m F + 2g t < l m F x fi Ns x m F ] . (10) 
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Eqs. dHllini) apply to the four ferromagnet/normal metal 
interfaces in a MJJ. We assume that the normal metal 
nodes have sizes smaller than the spin flip diffusion 
length, allowing its spin flip scattering to be neglected. 
Setting the total charge and spin-currents into each of 
the two normal metal nodes to zero yields eight equations 
for the eight unknown parameters which specify the spin- 
dependent chemical potentials of the two normal metal 
nodes. In the following two sections, we shall calculate 
the resistance for both the static regime and the dynamic 
regime from this set of equations. 



IV. STATIC MJJ RESISTANCE 

To apply magnetoelectric circuit theory to an MJJ 
with a time-independent free magnet, we note that m F — 
±z for source and drain electrodes respectively, set the 
source and drain chemical potentials to ±17/2 and assume 
that source and drain ferromagnetic electrodes are iden- 
tical. It follows from symmetry that spin-accumulation 
in the direction perpendicular to m and z vanishes in 
both normal metal nodes, reducing the number of free 
variables to six. We separate the normal metal chemi- 
cal potentials into accumulation (+) and transport (— ) 
contributions defined by 

Ma,± = Va,Nl ± fJ-a,N2, (H) 



electrodes fi z ,+ - 

2(g + 2G)fi c _ + g P fi Zj+ + 2GP^ TO ,_ = 2gU, 
2gpHc,- + (.9 + Grj)p, Zt+ = 2gpU, 
4GP/x c ,_ + (2G + W)Mm,- = 0, (12) 

where we have introduced the shorthand notations g = 
9 + 9 i P = (<7 — 9 )/<? and rf — 2g''^/g for the 
interfaces between normal metal and source and drain 
electrodes. The corresponding upper case G's and r\ re- 
fer to interfaces between the normal metals and the free 
magnetic layer. 

The charge current is most conveniently evaluated at 
either source or drain electrode with the result that 
I c = (eg/2h)((x c< - +PHz,+/2 — U). Although the static 
MJJ resistance Rs = U/(—eI c ) depends on the detailed 
values of the various interface conductance parameters, it 
tends to be close to the sum of the individual interface re- 
sistances and is not especially sensitive to source or drain 
electrode spin-polarization; in the limit in which all fcr- 
romagnets are fully polarized and identical, for example, 
the static MJJ resistance is 5(h/e 2 )/g. This result applies 
until the perpendicular spin current I z s = (Grj/Sir) fJ. z ,+ 
reaches its critical value J cr i t = K p . The corresponding 
critical charge current flowing through the MJJ is propor- 
tional to J C rit. For typical values p ~ 1, T) ~ 2, K p ~ 10 5 J 
(per m 3 ) and assuming that G = g, the citical charge cur- 
rent density is on the order of 1 x 10 -7 A/nm 2 for a free 
layer with a thickness around lnm. 



where a = c, m, z for charge, free magnet, and source 
spin directions, and Ni (i = 1, 2) denotes the i-th normal 
metal node. 

Simple equations for [i a ,± can be obtained by adding 
and subtracting the charge and spin conservation condi- 
tions for the two nodes. In this way it is easy to verify 
that for the static case /i c ,+ = /i m ,+ = Mz,- = 0, prop- 
erties that also can be established by symmetry consid- 
erations. The three remaining equations determine the 
chemical potential differences that drive charge and spin 
across the free magnetic layer, ^ c> _ and /i m ,-, and the 
accumulation of spin injected by the source and drain 



V. DYNAMIC MJJ RESISTANCE 

When l z s is substantially larger than / crit = K p a dy- 
namic steady state is approached in which 



(13) 



This precession of magnetization induced by spin current 
is the working principle of STT-oscillators. Because of 
their long spin relaxation times, the normal metal nodes 
do not follow this precession so that the components of 
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the spin-dependent part of the chemical potentials vanish 
along directions other than z. This leaves four unknown 
chemical potentials and <\> to be determined by Eq. (fT3"|) 
and the four circuit equations for conservation of I c and 
I z at each node. In the dynamic case the total spin- 
current I z s flowing out of the free magnet includes a spin- 
pumping contribution: 14 ' 17 



i; = ~G V (iM S , + -2t^). 



(14) 



The induced pumping current tends to counter the bias 
current, effectively increasing the resistance. Following 
the same strategy as in the static case we find that fJ. c ,+ = 
fj. z ,— — while the remaining unknowns satisfy 

2(5 + 2G)Mc-+fl»**,+ = 2 9 u > 
ZgPVc- + (.9 + Grj)fi Zt+ - 2Gr\Ti<\> = 2gpU, 

Gr)fj,z,+ - 2Gi]H = (15) 

The third of these equations relates the spin-current out 
of the free nanoparticle to the spin decay rate. The 
electrical properties of a MJJ are most interesting when 
the dimensionless coupling constant formed by compar- 
ing the right and left hand sides of Eq. ([15]) , 



A = 



AwaS 
HGr] ' 



(16) 



is small. Taking typical values^ of G ~ 10 (per 
nm 2 ), t) ~ 2 and of (4itM s )/j ~ 100 nm~ 3 yields 
A ~ 10ad[nm], where d is the film thickness, suggest- 
ing that small values are possible in thin magnetic films 
made from materials like permalloy that have small val- 
ues of a. Solving Eqs. (JT5J) for A = we find that the 
dynamic MJJ resistance 



_ h 4G + 2g( l-p 2 ) 
D e 2 4Gg(l 



P 2 ) 
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which diverges for p — > 1. When source and drain elec- 
trodes are strongly spin-polarized and A is small, the MJJ 
circuit should suffer a large increase in resistance at its 
critical current. 

This drastic resistance change can be qualitatively ex- 
plained as follows. In the static regime, the free layer 
magnetized in the sy-plane breaks spin rotation symme- 
try along z-axis. Exchange-field driven spin flips then 
provides a path for electrical conduction between oppo- 
sitely polarized electrodes. In the dynamic regime, pre- 
cession of the free layer restores spin rotation symmetry 
and the ^-component of the conduction electron spin is 
approximately conserved, increasing the electrical flow 
resistance. 

In Fig. [5J we plot the ratio £ = Rd/Rs, which char- 
acterizes the performance of a MJJ, as a function of the 
damping parameter a and the polarization p. £ increases 
as a decreases or p increases, and diverges in the limit 
a — > and p 1 . 
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FIG. 2: (color online). The ratio of dynamic and static re- 
sistance £ = Rd/Rs plotted as a function of damping pa- 
rameter a and polarization p. The parameters used here are 
g = G = 10 (per nm 2 ), P = 0.5, <q = 77' = 2, M s = 10 5 A/m, 
and d = lnm. This ratio £ diverges as a — > and p — > 1. 
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FIG. 3: (color online). The ratio between h<j>/e and the mag- 
netization dynamics induced emf in a MJJ plotted as a func- 
tion of polarization p for different values of damping param- 
eter q. These results were obtained for P = p, g = G = 10 
(per nm 2 ), r) = f]' = 2, M s = 10 5 A/m, and d = lnm. The 
straight line /3 — p applies for smooth magnetization textures. 



VI. DISCUSSION 



The influence of magnetization dynamics on electronic 
transport can be described by accounting for spin pump- 
ing, as in the analysis summarized above, or equivalcntly 
by accounting for a dynamics induced emf. For a slowly 
varying magnetic texture like a domain wall in transition 
metal ferromagnets, a transport electron follows the lo- 
cal magnetization direction adiabatically. The induced 
emf then given by pheb/e with <fi being the domain 
wall precession frequency and p the bulk transport cur- 
rent spin-polarization. For the MJJ considered here, spin 
transport is no longer adiabatic. In order to relate the 
two perspectives we assume that the influence of pre- 
cession in the MJJ can also be captured by imposing an 
extra emf fitvp/e that accounts for the difference between 
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static and dynamic circuit resistances: 

U _U -/3hj) 
Rd Rs 



(18) 



The effective value of ft then depends on microscopic MJJ 
parameters and can be evaluated by comparing Eq. (|T5|) 
with the results of a magnetoelectronic circuit analysis. 
Some typical values of ft are plotted in Fig.[3]as a function 
of p (= P) for different values of a. For the completely 
adiabatic case, we should expect ft = p which is a straight 
line in the figure. Although the ft of a MJJ generally de- 
viates from this straight line due to nonadiabatic effects, 
it approaches the straight line in the p — > and p — > 1 
limits. For p = 1, the same result can be obtained by 
making a unitary transformation into the rotating frame 
of rh and the result amounts to a relative shift hcj) be- 
tween the chemical potentials of the two electrodes. We 
find that the magnetization dynamics damping, parame- 
terized by a, slightly decreases the value of ft because the 
dissipation of energy through damping can be regarded 
as an extra resistance in the circuit. 

From the analysis above, the current drops once bias 
voltage U increases above a critical value determined by 
the in-plane anisotropy, i.e. when the precession of mag- 
netization begins. This signals a negative resistance re- 
gion of the FV curve similar to that for the Gunn diode. 19 
For MJJ with high performance (£ 3> 1), it could be used 



as a good fuse link, for example, to control the maximum 
current that can flow through a circuit. 

To summarize, we propose a ferromagnet/normal 
metal heterostructure which has properties analogous to 
those of a Josephson junction. We analyze its I-V charac- 
teristics and identify a static regime in which the magne- 
tization configuration is static with a low resistance and a 
dynamic regime in which the magnetization precesses and 
the resistance is high. Our analysis indicates that MJJ 
circuits which exhibit large and sudden relative change 
in resistance can be fabricated using modern spintronic 
materials with perpendicular anisotropy— and strongly 
polarized spin currents^. Besides having an academic 
interest, because of the insight that it provides in com- 
paring current-induced emf and spin-pumping notions in 
metal spintronics, the MJJ may be useful for applica- 
tions. 
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